Search for: All records

An organism becomes genetically mosaic through the accumulation of somatic mutations. Genetic mosaicism is a commonality of multicellular life and has been studied extensively in humans due to its associations with aging and diseases. In humans, somatic selection shapes the accumulation of somatic mutations, with strong signatures of positive somatic selection in cancer cell lineages. So far, evidence for somatic selection in plants has been inconsistent. The evolutionary implications of genetic mosaicism in humans and other animals are limited by early specification of germline cells, preventing transmission of somatic mutations to progeny. In contrast, many plant lineages reproduce asexually with clonal progeny derived from vegetative tissues. We describe the patterns and processes shaping somatic mutation accumulation within a single, 149-year-old historic sweet orange (Citrus sinensis) tree and within a clonal lineage of sweet orange. More than 12,000 somatic mutations were identified in the historic tree and 28,000 somatic mutations were identified across 199 clonally related sweet orange accessions. Both the spatial and genomic distributions of somatic mutations are non-random. The spatial patterns of somatic mutations across the historic tree depend on tree growth and development and their accumulation across the tree canopy recapitulates branching topology. Analysis of the genomic distribution of somatic mutations revealed that the subtelomeres, which are large arrays of ~180 bp repeats, are mutation hotspots. Finally, there was genomic evidence that somatic selection shapes the accumulation of somatic mutations both within the historic tree and also during clonal propagation. 

Interspecific hybridization is a common and effective strategy for producing disease resilient citrus cultivars, including those with tolerance to Huanglongbing (HLB) disease. Several HLB-tolerant cultivars have been developed through hybridization of mandarins (Citrus reticulata) with their wild relativePoncirus trifoliata. One such cultivar, ‘US-897’, exhibits robust tolerance to the bacteria causing HLB disease,Candidatus Liberibacter asiaticus(CLas). To explore the genetic architecture of the early transcriptional response toCandidatusLiberibacter asiaticus (CLas) infection in ‘US-897’, we performed transcriptomic analysis of the hybrid and its parents, ‘Cleopatra’ (C. reticulata) and ‘Flying Dragon’ (P. trifoliata). A haplotype-resolved genome for ‘US-897’ was generated using PacBio HiFi sequencing reads to support quantification of the expression of both theCitrus and Poncirusalleles. By profiling gene expression in this parent-offspring trio, we were able to determine the mode of inheritance for genes differentially expressed between parents (‘Cleopatra’ and ‘Flying Dragon’) and their interspecific hybrid (‘US-897’), with the majority genes exhibiting non-additive patterns of gene expression inheritance. Additionally, analysis of allele-specific expression in the hybrid ‘US-897’ revealed the contribution of cis- versus trans-acting regulatory variants on genes with additive and non-additive modes of inheritance. A strong correlation between differential expression between parents and allele-specific expression in ‘US-897’ suggests that cis-regulatory variation is a significant source of expression divergence between species. Finally, genes responsive to infection withCLas were identified to explore how gene regulation associated with tolerance to HLB was rewired betweenCitrusand its relativePoncirus. 

Abstract BackgroundGenetic and epigenetic perturbation of cis-regulatory sequences can shift patterns of gene expression and result in novel phenotypes. Phased genome assemblies now enable the local dissection of linkages between cis-regulatory sequences, including their epigenetic state, and allele-specific gene expression to further characterize gene regulation and resulting phenotypes in heterozygous genomes. ResultsWe assembled a locally phased genome for a mandarin hybrid named ‘Fairchild’ to explore the molecular signatures of allele-specific gene expression. With local genome phasing, genes with allele-specific expression were paired with haplotype-specific chromatin states, including levels of chromatin accessibility, histone modifications, and DNA methylation. We found that 30% of variation in allele-specific expression could be attributed to haplotype associated factors, with allelic levels of chromatin accessibility and three histone modifications in gene bodies having the most influence. Structural variants in promoter regions were also associated with allele-specific expression, including specific enrichments of hAT and MULE-MuDR DNA transposon sequences. Integration of haplotype-resolved genetic and epigenetic landscapes with high-throughput phenotypic analysis of fruit traits in a panel of 154 accessions with mandarin and pummelo ancestry revealed that trait-associated variants were enriched in regions of open chromatin. Mining of trait-associated variants uncovered a Gypsy retrotransposon insertion in a gene that regulates potassium transport and may contribute to the reduction in fruit size that is observed in mandarins. Conclusions​​Using a locally phased assembly of a heterozygous cultivar of citrus, we dissected the interplay between genetic variants and molecular phenotypes to reveal cis-regulatory sequences with potential functional effects on phenotypes relevant for genetic improvement. 

Abstract Genomic clusters of immune genes, including those encoding nucleotide-binding leucine-rich repeat (NLR) proteins, are a model for exploring the dynamics of genomic regions in flux. Rapid sequence evolution of immune genes, including NLRs, and variation in their gene content, may enable long-lived plants, which lack adaptive immune systems, to keep pace with the fast evolution of pathogens. To explore the patterns and processes shaping the evolution of NLR gene content in a genus of long-lived tree species, we unified the annotation of NLR genes across 11 accessions (or 15 haplotypes) from the genusCitrusand its relatives, including three new diploid genome assemblies. A majority of NLRs were arranged in genomic clusters composed of paralogous genes, typically from a single gene family. Even larger clusters, with 10 or more NLRs, were limited to genes derived from one or few gene families. These patterns suggested that genomic clustering of NLRs arose through local expansion of phylogenetically related NLRs, but the mechanistic processes driving these patterns are not clear. Local gene duplication can be mediated by multiple processes, including transposon-mediated gene capture and subsequent proliferation, and non-allelic repair of double stranded breaks, including unequal recombination. Examples of retrotransposon-mediated duplication of NLRs were identified, but these were not sufficient to explain massive regional expansions. Signatures of unequal recombination are challenging to identify. Focusing on recent lineage-specific sequence duplications, at least one case of unequal recombination was identified, supporting a role for unequal recombination in shaping genomic variation in these regions. 

Abstract Gene body methylation (gbM) is an epigenetic mark where gene exons are methylated in the CG context only, as opposed to CHG and CHH contexts (where H stands for A, C, or T). CG methylation is transmitted transgenerationally in plants, opening the possibility that gbM may be shaped by adaptation. This presupposes, however, that gbM has a function that affects phenotype, which has been a topic of debate in the literature. Here, we review our current knowledge of gbM in plants. We start by presenting the well-elucidated mechanisms of plant gbM establishment and maintenance. We then review more controversial topics: the evolution of gbM and the potential selective pressures that act on it. Finally, we discuss the potential functions of gbM that may affect organismal phenotypes: gene expression stabilization and upregulation, inhibition of aberrant transcription (reverse and internal), prevention of aberrant intron retention, and protection against TE insertions. To bolster the review of these topics, we include novel analyses to assess the effect of gbM on transcripts. Overall, a growing body of literature finds that gbM correlates with levels and patterns of gene expression. It is not clear, however, if this is a causal relationship. Altogether, functional work suggests that the effects of gbM, if any, must be relatively small, but there is nonetheless evidence that it is shaped by natural selection. We conclude by discussing the potential adaptive character of gbM and its implications for an updated view of the mechanisms of adaptation in plants. 

Abstract Methylated CHH (mCHH) islands are peaks of CHH methylation that occur primarily upstream to genes. These regions are actively targeted by the methylation machinery, occur at boundaries between heterochromatin and euchromatin, and tend to be near highly expressed genes. Here we took an evolutionary perspective by studying upstream mCHH islands across a sample of eight grass species. Using a statistical approach to define mCHH islands as regions that differ from genome-wide background CHH methylation levels, we demonstrated that mCHH islands are common and associate with 39% of genes, on average. We hypothesized that islands should be more frequent in genomes of large size, because they have more heterochromatin and hence more need for defined boundaries. We found, however, that smaller genomes tended to have a higher proportion of genes associated with 5′ mCHH islands. Consistent with previous work suggesting that islands reflect the silencing of the edge of transposable elements (TEs), genes with nearby TEs were more likely to have mCHH islands. However, the presence of mCHH islands was not a function solely of TEs, both because the underlying sequences of islands were often not homologous to TEs and because genic properties also predicted the presence of 5′ mCHH islands. These genic properties included length and gene-body methylation (gbM); in fact, in three of eight species, the absence of gbM was a stronger predictor of a 5′ mCHH island than TE proximity. In contrast, gene expression level was a positive but weak predictor of the presence of an island. Finally, we assessed whether mCHH islands were evolutionarily conserved by focusing on a set of 2,720 orthologs across the eight species. They were generally not conserved across evolutionary time. Overall, our data establish additional genic properties that are associated with mCHH islands and suggest that they are not just a consequence of the TE silencing machinery. 

Abstract In plants, mammals and insects, some genes are methylated in the CG dinucleotide context, a phenomenon called gene body methylation (gbM). It has been controversial whether this phenomenon has any functional role. Here, we took advantage of the availability of 876 leaf methylomes in Arabidopsis thaliana to characterize the population frequency of methylation at the gene level and to estimate the site-frequency spectrum of allelic states. Using a population genetics model specifically designed for epigenetic data, we found that genes with ancestral gbM are under significant selection to remain methylated. Conversely, ancestrally unmethylated genes were under selection to remain unmethylated. Repeating the analyses at the level of individual cytosines confirmed these results. Estimated selection coefficients were small, on the order of 4 Nes = 1.4, which is similar to the magnitude of selection acting on codon usage. We also estimated that A. thaliana is losing gbM threefold more rapidly than gaining it, which could be due to a recent reduction in the efficacy of selection after a switch to selfing. Finally, we investigated the potential function of gbM through its link with gene expression. Across genes with polymorphic methylation states, the expression of gene body methylated alleles was consistently and significantly higher than unmethylated alleles. Although it is difficult to disentangle genetic from epigenetic effects, our work suggests that gbM has a small but measurable effect on fitness, perhaps due to its association to a phenotype-like gene expression. 

Plant genomes interact when genetically distinct individuals join, or are joined, together. Individuals can fuse in three contexts: artificial grafts, natural grafts, and host–parasite interactions. Artificial grafts have been studied for decades and are important platforms for studying the movement of RNA, DNA, and protein. Yet several mysteries about artificial grafts remain, including the factors that contribute to graft incompatibility, the prevalence of genetic and epigenetic modifications caused by exchanges between graft partners, and the long-term effects of these modifications on phenotype. Host–parasite interactions also lead to the exchange of materials, and RNA exchange actively contributes to an ongoing arms race between parasite virulence and host resistance. Little is known about natural grafts except that they can be frequent and may provide opportunities for evolutionary innovation through genome exchange. In this review, we survey our current understanding about these three mechanisms of contact, the genomic interactions that result, and the potential evolutionary implications.